Non-adiabatic Fast Control of Mixed States based on Lewis-Riesenfeld Invariant 
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Abstract 

We apply the inversely-engineered control method based on Lewis-Riesenfeld invariants to control mixed states of a two-level 
quantum system. We show that the inversely-engineered control passages of mixed states - and pure states as special cases - can 
be made significantly faster than the conventional adiabatic control passages, which renders the method applicable to quantum 
computation. We devise a new type of inversely-engineered control passages, to be coined the antedated control passages, 
which further speed up the control significantly. We also demonstrate that by carefully tuning the control parameters, the 
inversely-engineered control passages can be optimized in terms of speed and energy cost. 

PACS numbers: 37.10.Dc, 32.80.Qk, 42.50.-p, 03.65.Ca, 03.65.Vf, 03.65.Yz, 03.67.-a, 03.65.Ta 
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I. INTRODUCTION 



Adiabatic control has been one of the main meth- 
ods of quantum control, quantum information processing, 
and quantum computation, because it is able to evolve 
a quantum system against parameter fluctuations. In 
particular, Farhi et al. proposed the adiabatic quantum 
computation a new model of quantum computation 
based on adiabatic control, which was then extended to 
quantum open systems 0, Q. Note, however, that adia- 
batic control gains its robustness on the cost of losing the 
speed of control, which seriously limits its applicability 
to cases where fast control - in particular faster than the 
decoherence of the system - is of great concern, such as 
quantum information processing and quantum comput- 
ing UM- 

In an adiabatic control of a quantum system, the sys- 
tem remains in the instantaneous ground state of its 
time-dependent Hamiltonian during the entire evolution. 
The control parameters in the Hamiltonian are care- 
fully designed such that the adiabaticity condition al- 
ways holds, which usually results in very long execution 
time. While other ways of speeding up an adiabatic con- 
trol arc available in the adiabatic regime Chen et 
al. put forward a new method, the inversely-engineered 
control (lEC) based on Lewis-Riesenfeld invariants [lol - 
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Il5l | . Lewis-Riesenfeld invariants (LRI) were discovered by 
Lewis and Riesenfeld in 1969 to solve the time-dependent 
Schrodinger equations [l^. In lEC, the Hamiltonian, 
which controls the system evolution, is determined by a 
corresponding LRI, which is constructed according to the 
desired initial and final states of the system. The lEC of 
a system is not adiabatic in general, although the asso- 
ciated Hamiltonian of the system may happen to satisfy 
the adiabaticity condition. Chen and his collaborators 
have shown that it was plausible that lEC offered short- 
cuts to adiabatic control of the pure states of a two-level 
quantum system without heating the system [Til [l^ . 

In this paper, we study lEC in more details by apply- 
ing it to mixed states of a two-level quantum system. We 
consider mixed states because they are more general than 
pure states and play the key role in NMR (Nuclear Mag- 
netic Resonance) quantum computation. In fact, how 
to control mixed states efficiently and rapidly is a ma- 
jor problem in NMR computation. We show that the 
lEC does give rise to fast, nonadiabatic control passages 
of a two-level system, compared with the viable adia- 
batic control passages of the system. Moreover, the lower 
bound of the control time of lEC is only constrained phys- 
ically but not mathematically. In particular, lEC bears 
an optimization regarding the time and energy cost of 
the control. 

The paper is organized as follows. We briefly review 
the LRI and lEC in the next section. Section HI is the 
main part of this paper, where the lEC of a mixed state 
is investigated. Section IV is devoted to conclusions and 
outlook. Throughout the paper, h is set to one. 
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II. LEWIS-RIESENFELD INVARIANT AND 
INVARIANT-BASED INVERSE ENGINEERING 

A Lewis- Riesenfeld Invariant I(t), defined as a dynami- 
cal invariant of a quantum system witli a time-dependent 
Hamiltonian Hit), is a Hermitian operator defined on tlie 
Hilbcrt space of tfie system, wliicli satisfies tlie following 
defining equation [l6j . 



.dm 

dt 



[H{t)J{t)]=Q. 



(1) 



0, 



The dynamical invariancc of I{t) implies d(/(t))/dt 
which endows I{t) with a spectral decomposition 



where A„ is real and time- independent. A LRI of a sys- 
tem enables the direct integration of the time-dependent 
Schrodinger equation of the system, giving rise to the 
general solution in terms of the instantaneous eigenstates 
of /(t) [H, i.e.. 



(2) 



to a group of ordinary differential equations of the param- 
eters in I{t) and those in H{t). 4) Plug in the specified 
initial and final conditions on the parameters in I{t). 5) 
Find physically and mathematically reasonable interpo- 
lations of the parameters in /(i), such that they comply 
with the initial and final conditions. 6) Finally, determine 
the parameters in H{t) by solving the group of differen- 
tial equations. The Hamiltonian H{t) determined this 
way is the desired one. This is indeed the strategy taken 
in [111 [T^ . Sarandy et al. adopts another strategy based 
on Lie algebras [l7| . 



III. INVERSELY-ENGINEERED CONTROL OF 
THE MIXED STATES OF TWO-LEVEL ATOM 

From now on, we focus on a two-level quantum system 
to study lEC in details. Explicitly, we have a spin-1/2 
nucleus under a strong magnetic field and an oscillating 
transverse radio frequency (RF) fields in our mind. This 
system is realized in NMR for example. In Sections IIIIl l 
to 3.3, we shall explicitly follow the steps of our strategy 
described in the end of the last section. 



where a„ is the Lewis-Riesenfeld phase 

«n = j\<i^n{t')\ {i^, - H{t')^ \Mmt'. (3) 

In an adiabatic control of a system with a Hamiltonian 
Hq, one usually applies external fields on the system, such 
that the total Hamiltonian becomes H{t), and drives the 
system from an initial Hamiltonian H{Q) = Hq to the 
final one, H{tf), by keeping the system in the instan- 
taneous ground state of H{t). The corresponding final 
state is either an eigenstate or a superposition of the 
eigenstates of H{tf), depending on the purpose of the 
control. We know that associated to this Hamiltonian is 
a LRI I{t), which satisfies Eq. ([l]). The converse is also 
true. That is, there exists a H{t) that corresponds to a 
given I(t), which indicates a way around the adiabatic- 
ity constraint. That is, rather than evolving a system 
adiabatically, one first constructs a LRI I{t), which then 
yields a Hamiltonian H{t) by Eq. ([T|), that drives the sys- 
tem in an instantaneous eigenstate \(j>n{t)) of I{t). States 
\(j)n{0)) and \(j)n{tf)) must coincide with the designated 
initial and final states of the original adiabatic control. 
At any intermediate time t, however, |(/)„(i)) is not an 
eigenstate of H(t) in general. This summarizes the idea 
of lEC. 

Here is a strategy to implement this idea on any sys- 
tem. 1) Write down the total Hamiltonian H{t) of the 
system and the control field, which reduces to at t = 
and t = tf, but keep all the time-dependent parameters 
undetermined. 2) Use the most convenient LRI I{t) of 
H{t), leaving the parameters undetermined either. 3) 
Plug both H{t) and /(<) into Eq. which should lead 



A. The Lewis-Riesenfeld Invariant 

The typical two-level NMR Hamiltonian under an ex- 
ternal RF field reads 



H{t) = - 



^gi[(w-|-tjRF)t-<^]l 



^C-i[("+'^RF)t-¥'] 



J 

(4) 

where flu, wrf, wq, w, and ip are respectively the Rabi 
frequency, the frequency of the RF field, the Zeeman en- 
ergy of the system, the frequency of the rotating frame, 
and a phase. It is customary to configure the system such 
that Lu+ujxp ^ oj— wrf, which leads to the rotating- wave 
approximation, in which the highly oscillating terms due 
to cj -f wrf ill Eq. (HI) are averaged to vanish. Moreover, 
one can put lu = wrf • 

For the purpose of control, however, we can modulate 
the Rabi frequency, the w, and the phase (p with time. 
Therefore, by defining the detuning frequency A{t) = 
Lo{t) — wq, the Hamiltonian of our concern hereafter is 



H = - 



nRit)e 
-A(t) 



iip{t) 



(5) 



This Hamiltonian has eigenvalues E± ~ zthn/2, where 
H. = y/ A'^ + ilj^ is the generalized Rabi frequency. At 
t = and t ~ tf when the RF field is turned off, i?(0) = 
~H{tf) ^ Ho = diag(A(0),-A(0)). In the sequel, we 
assume ip{t) = 0, which further simplifies the H{t) and 
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renders the adiabaticity condition of the system 
UrA - iln A 
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(6) 



and the instantaneous eigenstates of the Hamiltonian 



|£;+(0)=cos( - J|0)+sin^-)|l), 



|i?_(t))=sin(-J|0)-cos^-)|l), 



(7) 



where 9 = arccos(A/r2) is the mixing angle. |0) = (1, 0)-'" 
and |1) = (0, 1)^ are the eigenstates of the Pauh matrix 
az- The Hamihonian H{t) in Eq. (O may serve to control 
the system adiabatically if it meets the adiabaticity con- 
dition in Eq. (jS]). Note that none of the time-dependent 
parameters in H{t) is fixed. This concludes the step 1) 
of our strategy. 

A reasonable LRI I{t) associated with H{t) in Eq. ([5]) 
must be employed to complete the step 2) of our strategy. 
We adopt the general form employed in [iH, i.e., 



cos 7 
e~^^ sin 7 



e*^ sin 7 
' cos 7 



(8) 



where 7 and /? are time-dependent parameters that in 
turn determines ^R{t) and A(t) in the Hamiltonian, 
while flo is an auxiliary constant. 

To follow step 3), we substitute Eqs. ^ and ([5]) into 
Eq. ([T]), which immediately leads to two coupled differ- 
ential equations 



7 



smp 

A = rtfj cot 7 cos (3 — $■ 



(9) 



In lEC, wc need first to find 7 and /3 without solving 
Eq. © and then insert them into Eq. © to obtain Hr 
and A. We now proceed to do this to control the mixed 
states of our two-level system. 



B. Inversely-Engineered Control of Mixed States 

To complete the step 4) , we first find the instantaneous 
eigenstates of I{t), as the two- level system evolves in a 
mixed state in the basis of the instantaneous eigenstates 
of I{t). From Eq. the instantaneous eigenstates of 
I{t) are obtained as 



|</>+(t))=cos(^)e^^|0)+sin(^)|l), 
|0_(t))=sinQ)|O)-cosQ)e-'^|l). 



(10) 



One can verify that in this case, the Lewis-Riesenfeld 
phase in Eq. ([3]) becomes 
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a±{t) = T:, / n{t')dt', 



(11) 



where r2(t) = {A + $) cos 7 + $ + sin 7 cos /?. Hence, 
the instantaneous mixed state is 



(12) 



1 1 + (P+ -P-)cos7 {p+-p- 



2 \ip+- p- 



sm7 

sin 7 1 — (p+ — p„) COS7 



where J2i=±Pi ^ 1- Note that the density matrix is 
diagonal in the basis of the eigenstates \(j)±) of I{t) but 
not in the basis \E±{t)) in general. Since we would like 
to depict the trajectories of the evolution of the mixed 
states in the Bloch ball for illustrative purposes, we cast 
the density matrix (t) in the basis of Pauli matrices as 

p'{t) =\h + ^i-y^[sin7(t)cos/?(t)a, 

— sin7(i) sin P{t)(jy + cos^{t)(Tz\^ (13) 

where I2 is the 2x2 unit matrix. The 
corresponding Bloch vector is {pj^ — p_) x 
(sin7(t) cos/3(t), — sin7(i) sin/3(t), cos7(t)). Clearly, if 
cither p+ or p_ is zero, {i) reduces to a pure state and 
the Bloch vector sits on the Bloch sphere. 

We shall consider the control passages of p^ that cause 
a complete population inversion in the two levels of the 
atom. Seen in Eq. (fTO)) . the population is completely 
controlled by 7 because 13 is merely a phase. There- 
fore, in order to achieve the population inversion, 7 
must satisfy certain initial and final conditions, which, 
following from Eq. (fTO)) . should be 7(0) = titt and 
7(t/) = mTT, where n and m are respectively even and 
odd integers, or vice versa. Without loss of general- 
ity, we assume 7(0) = tt and 7(t/) = 0, such that 
p^(0) = diag((l -P++ p_)/2, {l+p+~ P-)/2)) and 
p'{tf ) = diag((l +P+ -p_)/2, (1 -p+ +p_)/2)), which 
are indeed the corresponding mixed states in terms of 
the two levels of Hq, manifesting the population inver- 
sion. Although the initial and final density matrices are 
diagonal, the intermediate ones are generally not. We 
can constrain 7 further because the RF field and hence 
the Rabi frequency Q,r is turned on at t = and off at 
t = tf, which sets 7(0) = 7(i/) = by Eq. Summa- 
rized in below are the necessary conditions for ^{t) 



7(0) = ^, 7(i/) = 0, 
7(0) = 0, ^{tf) = 0. 



(14) 



We now examine /3. As previously remarked, (3 plays 
no role in the evolution of the population number. How- 
ever, it enters the Bloch vector corresponding to p^ and 
thus affects the evolution trajectories of the states. Fur- 
thermore, Eq. © shows that (3 also influences the profiles 
of and A. That is, inappropriate j3 may make VIr and 
A unphysical — ruins the rotating wave approximation, 
for example — and may cause the control impractical. 
We fix two most obvious conditions for (3 for the moment 
but leave other conditions for the subsequent sections 
because they are case-dependent. Since 7(0) = tt and 



3 



7(t/) = 0, A can diverge at t = and t = tf because 
of the cot 7 in Eq. ^ unless this divergence is compen- 
sated by vanishing cos /?. This cahs for the following two 
conditions. 

m = I3{tf) = -tt/2. (15) 
They are negative in order that lim f2/j ~ 0"*", which is 

t~¥Q,tf 

consistent with the positivity of 51 due to the assump- 
tion that if = 0. 

It is time to proceed to step 5) to find parameters ^(t) 
and (3(t) that satisfy the above conditions and certain 
extra conditions to be set. To this end, Ansatz about 
the functional forms of 7 and /3 should be made. The 
simplest one is that both 7 and /3 are polynomials of fi- 
nite order in t. We shall take this ansatz for our study in 
all following sections but one, where it will be explicitly 
stated. It turns out that the power of lEC in speed- 
ing up the control of the system varies in different cases. 
Throughout the rest of the paper, we take p+ = 0.2 and 
P- = 0.8 for dcfiniteness. Our analysis, however, applies 
to any states including pure states. 
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FIG. 1: (Color online.) lEC (dashed lines) and adiabatic 
control (solid) passages of the mixed state in Eq. (|12|l . tf — 1. 



final values of the detuning, such that H{tf) = —H{0). 
Plugging 7 and /3 so obtained into Eq. (|T^ gives the cor- 
responding p^{t) with t G [0,i/]. The solid curves in Fig. 
[T]show the control passages p{i{t) and P22{t)- Step 5) is 
now completed. 

We remark that in Fig. [T] and in all subsequent fig- 
ures, the final time tf can have any positive value, such 
as microseconds or nanoseconds. In other words, math- 
ematically, tf can be arbitrarily close to zero. Physi- 
cally, of course, exceedingly small tf may set the en- 
ergy/frequency scale of fin and A too high to be practical 
or realizable. Another physical lower bound oi tf roots 
in the uncertainty principle [l8l | . 
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FIG. 2: (Color online.) Bloch vector trajectories of the mixed 
state under the lEC (outside the xz-plane) and the adiabatic 
control (lying in the xz-plane), corresponding to Fig[l] Three 
dots represent respectively the sphere center and where the 
lEC and adiabatic trajectories intersect the ly-plane. The 
solid straight lines connecting the three dots help with the 
visualization of the spatial relations of the two trajectories. 
The shaded disk in the sphere is the xz-plane. 



C. Third Order 7 

If we consider the conditions on 7 in Eq. ((Ti)) exclu- 
sively, then a third-order polynomial j{t) — X]i=o "^^^^ 
a good candidate for the functional form of 7 that inter- 
polates these conditions. Solving for the coefficients a^'s 
by substituting 7 into Eq. shows that 7(i) monoton- 
ically decreases from tt to with t G [0, t/]. This negativ- 
ity of 7(i) indicates that sin/3 must be negative for flfj 
to remain positive by Eq. ([9]). As such, we assume /3(t) 
is also third-order in t and accordingly add to Eq. 
two more conditions on /3, i.e., ${0) = 3n/{2tf) and 
l3{tf) = — 37r/(2i/)[2l|, which simply fix the initial and 



The parameters Clj^ and A calculated from 7 and /3, 
corresponding to the control passages in Fig. [1] happen 
to satisfy the adiabaticity condition Eq. ([5]). Hence, one 
can directly perform an adiabatic control of the mixed 
state without recourse to LRI. The adiabatically evolving 
density matrix reads 

p^'{t) = Y,p,\E,it)){m)\ (16) 

= + ^ [sin +cos0az], 

where the second line is the Bloch vector decomposition, 
and 9 = 9{t) is the mixing angle. The matrix elements 
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Pif{t) and P22{t) of the adiabatic control passages are 
depicted as solid lines in Fig. [TJ Note that they almost 
coincide with the lEC passages. In this case, the lEC is 
as slow as the adiabatic control. Nonetheless, the tra- 
jectories of the mixed state evolving under the lEC and 
adiabatic control respectively are rather different in the 
Bloch ball as visualized in Fig. [2] It is clear from Eqs. 

and that the adiabatic Bloch vector is confined 
within the a;2;-plane but the lEC one is not. 

In theory, the Hamiltonian inversely-engineered from 
a LRI may violate the adiabaticity condition Eq. ([6]). In 
such a scenario, if one still insisted on performing an 
adiabatic control with this particular Hamiltonian, one 
would have to accomplish it in literally an infinite amount 
of time. Consequently, when the inversely-engineered 
Hamiltonian disallows an adiabatic control in an accept- 
able time duration tf, lEC outmatches adiabatic control 
significantly because its control time is bounded below 
only by the physical realizability of the control parame- 
ters and the uncertainty principle. This motivates us to 
investigate how to accelerate lEC more and optimize it in 
the nonadiabatic regime. Recall that /3 does not directly 
affect the evolution of the population. Thus, we should 
first manipulate 7 and then /3 if necessary. 



D. Fourth Order 7 

We did not have much freedom in manipulating 7 in 
the case where it was a third-order polynomial in t. It 
follows naturally to assume fourth-order 7, namely 7(i) = 
Ylt=Q^i^^ ■ coefficients in 7, one must 

add an extra condition to Eq. We impose conditions 
on the value of ^{tf /2) as 

(i) 7(i//2) - 27r/5, (ii) 7(i//2) = 27r/6. (17) 

Let us first consider the condition (i) in Eq. (|T7)) . The 
fourth-order 7. with this new condition, again decreases 
monotonically from tt to as i changes from to tf. 
The regular behavior of the new 7 enables us to reuse 
the third-order /? obtained in the previous section. The 
corresponding lEC gives rise to control passages as shown 
by the dot-dashed curves in Fig. 121 which are obviously 
faster than the dashed ones because the fidelity of the 
population inversion in the former approaches one faster. 
The latter ones are in fact the lEC passages in Fig. [T] 
Therefore, this new fourth-order 7 expedites the lEC. 

In fact, further acceleration can be achieved by us- 
ing the condition (ii) in Eq. ()17p instead, which assigns 
7(t//2) a value smaller than that in the condition (i) in 
Eq. ([T7]). This results in the control passages shown as 
the dashed curves in Fig. [31 Figure [H depicts the and 
A corresponding to the solutions in Fig. [3l One might 
wish to further speed up the control by tuning ^{tf/2) 
down continuously to some finite and positive value. 
Unfortunately, this is infeasible because the condition 
7(t//2) = 27r/6 is already near the limit, 27r/6. 40175, 
below which the fourth-order 7 that meets the condition 
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FIG. 3: (Color online.) lEC control passages for different 
intermediate conditions in the legend on 7 at tf/2. 



behaves oddly in the sense that 7(t) has a sharp drop 
into the negative domain at time near tf. If /3 remains 
third-order, singular behavior occurs in A. The worse is 
that one cannot design a different (3 that is able to di- 
minish the singularity in A without introducing singular 
behavior in CIr. 




FIG. 4: (Color online.) Parameters (curves completely 
above the t — axis) and A (curves intersecting the t = 
axis) associated with the control passages in Fig. [3] respec- 
tively, with the same plot style. 

Nevertheless, there exist cases of lEC where certain 
singular behaviors in flu and/or A caused by 7 can be 
remedied by /3 if /3 has some special profiles. When this 
happens, the lEC time significantly diminishes if another 
technique, to be introduced shortly, is employed. We 
elucidate this by an example. Suppose we take 7(t//2) = 
and add this condition to Eq. (fT4l) to impose a new set 
of conditions on 7 as 



(i) 7(0) = TT, (ii) Jitf) = 0, 

(iii) 7(0) = 0, (iv) jitf) = 0, 
(v) ^{tf/2) = 0. 



(18) 
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Conditions in Eq. ([T5)) yield a fourth-order solution of 7, 
as plotted in Fig. [S] (dashed curve). 
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FIG. 5: (Color online.) Parameters 7 (dashed) and /3 (solid) 
as solutions respectively to Eqs. psp and p9p . We set tf — 1. 

As seen in Fig. [SJ 7(i) reaches zero at tf/2, and 7(i) 
changes sign at llt//16. The former makes cot 7 blow 
up at tf/2, while the latter may flip the sign of flu by 
Eq. Nonetheless, these issues can be made harmless 
by designing /? that satisfies the following conditions. 

(i) m = -n/2, (ii) I3{tf ) - tt/2, 

(iii) Pitf/2) - -7r/2, (iv) /3(1K//16) = (19) 

(v) /3(0) = -/3(i/)=V2i/. 

Condition (iii) in Eq. ([T5| keeps the term cot 7 cos (3 in 
Eq. dni) finite at tf/2 despite 7(i/) = 0. Notice from 
Eq. that limt^o 7 = 0^ and limj_i.j^72 7 = 0+. Con- 
ditions (i) and (ii) in Eq. ([T^ make 7/sin/3, i.e., flji, 
approach zero from above in these cases. Because /3 
changes sign with 7 at llt^/16, ^Ir remains nonnega- 
tive over all time. Moreover, condition (iv) in Eq. 
and that 7(llt//16) = ensure £7^ be finite at lli//16. 
Two more conditions should present themselves to fix the 
detuning at t = and t = tf, which for this example we 
choose to be Eq. ([TO]) (v). This particular choice seems 
to have no a priori reason over all possible choices. In- 
deed, as to be seen in Section IIIII5. the value of ${0) 
plays a key role in optimizing the energy cost of lEC, 
and Eq. ([T^ (v) is not the optimal choice. 

The sohd curve in Fig. [S] depicts the solution of /3 under 
the conditions Eq. (|T9| . which demonstrates the expla- 
nations just given. The 7 and /? in Fig. [5] produce flji 
and A as the solid curves in Fig. U and the corresponding 
control passages (solid) in Fig. |31 Interestingly, the solid 
control passages in Fig. [3] appears to attain population 
inversion with fidelity 1 at t ~ tf/2, much earlier than 
all others in the figure, notwithstanding the fidelity drops 
after t = tf/2 and re-approaches one ed, tf. 

Does this surprising phenomenon truly give the pop- 
ulation inversion as that to be achieved at tfl The 



answer is "Yes", provided that anther technique is em- 
ployed, which is switching the total Hamiltonian H{t) at 
t ~ tf/2 sharply to some constant one in which the pop- 
ulation inversion is well defined. Population inversion of 
a two-level system is well-defined only by the Hamilto- 
nian of the system without external fields. As such, the 
basis of the initial and final states of the control must 
be the eigenstates of Hq. Since sin[/3(t/)] is finite, condi- 
tions (iii) and (iv) in Eq. ([T5)) sets ri_R(0) and V,fi{tf) to 
zero, so that [7J(0),/(0)] = [H {t f) , I {t f)] = 0. Accord- 
ingly, each eigenstate of /(O) and that of I{tf) are also an 
eigenstate of H{Q) and one of Hitf), respectively. As we 
began the control with the state p^(0) = diag((l — p+ + 
P-)/2, (1 + p+ — p„)/2)), if we follow the solid control 
passage in Fig. [3] all the way to tf, we end up with the 
state p'{tf) = diag((l +P+ -p_)/2,(l -p+ +p_)/2)), 
which is the very population-inverted mixed state in the 
two- level basis of i/o- 

If, on the other hand, the evolution halts at t = tf/2 
along the solid control passage in Fig. [3l the mixed state 
in the eigenbasis of I{tf/2) is p^{tf/2) = diag((l +p+ — 
P-)/2, {\—p++p-)/2)), the same as p^{tf). Nonetheless, 
we only impose 7(t//2) = in Eq. (|T5)) but not ^{tf/2) = 
0. Thus, [H{tf/2),I{tf/2)] ^ 0, and p\tf/2) is not the 
right population-inverted state in the eigenbasis of the 
Hamiltonian. 

Fortunately, this issue has a resolution. In many ap- 
plications of quantum control, such as in NMR, the to- 
tal Hamiltonian H{t) is actually piecewise constant or 
piece wise continuous. Therefore, at tf/2, one can stop 
the control passage of the system and instantaneously 
switch the total Hamiltonian to a constant diagonal 
Hamiltonian in order that it commutes with I{tf/2), 
which is also diagonal. To do so, we simply need to 
turn off the RF field right at t = tf/2, i.e., setting 
r2i{(t > tf/2) = 0, then the Hamiltonian becomes 
H{tf/2) = diag(A(t//2), -A(i//2)), which shares the 
basis with p^{tf/2); hence, a well-defined population in- 
version completes. Note that A{tf/2) ^ A(t/) in gen- 
eral but this does not affect the definition of population 
inversion unless A{tf/2) = 0. Therefore, the popula- 
tion inversion at tf/2 along the solid control passage in 
Fig. [3] should be authentic, which exhibits a tremendous 
decrement in the control time. We christen the control 
passages such as the solid curves in Fig. [31 owing to con- 
ditions as in Eq. (|18p and alike, the antedated control 
passages to distinguish them from the usual ones that 
end at </. 

One may wonder why the antedated control passages 
are meaningful, as one can simply set the control time 
of the usual control passages to be some t'^ < tf, say, 
t'j: = tf/2, given the aforementioned fact that tf can 
be anything mathematically. Here is our reasoning. On 
the one hand, for whatever tf chosen for the usual con- 
trol passages, the antedated ones can always complete 
the control much earlier. On the other hand, the energy 
cost due to the antedated control passages can be made 
smaller than that due to the usual ones. Actually, all 
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the usual lEC passages in Fig. [3] have the same energy 

cost, namely /J^ nji{t)dt ~ 3.1482. One may keep this 
number in mind to compare it with the energy cost of 
the antedated control passages introduced in the next 
section. 

The antedated control passages speed up the control 
at a price of weakened stability. As seen in Fig. [31 the 
population inversion achieved by the antedated control 
passages (solid) appears at the peak and valley of the 
passages; hence, as opposed to the case of usual control 
passages, the evolution of the state must halt sharply at 
tf /2 because somewhat off that specific temporal position 
would reduce the fidelity. This price certainly depends 
on how much tolerance in compromising the fidelity is 
practical. 

We remark that although switching between Hamilto- 
nians takes some amount of time, turning on or off a RF 
field in NMR costs neg ligible time compared with the 
control time of lEC [13 • 

In the next subsection, we will look for more antedated 
control passages that are even faster and study their op- 
timization in terms of control time and energy cost. 

E. Optimizing lEC 

Condition (v) in Eq. urges one to modify it by 
merely shifting the antedated time ta, at which 7 hits 
zero, to attain faster antedated control passages. The 
time ta can be even earlier than tf/2. Instead of the 
condition (v) in Eq. (jlSp . we consider three different in- 
termediate conditions: 

(v')7(i//3)=0, (v")7(2t//7) = 0, (v'") 7(2t//7.6) = 0. 

(20) 




FIG. 6: (Color online.) Antedated control passages due to 
fourth-order 7 and fifth-order /3. The antedated time of each 
passage is shown in the legend. 

We take (3 fifth-order in t and require it to satisfy con- 
ditions adapted to each case in Eq. ([^0)1 separately, sim- 



ilarly to what we did to obtain Fig. [51 Basically, con- 
ditions (iii) and (iv) in Eq. (|T^ should be modified in 
accordance with the three cases in Eq. (|20|) . Nonetheless 
we choose not to show the new conditions on /? here, as 
they are irrelevant for our purpose at the moment. The 
antedated control passages obtained are shown in Fig. [51 
which demonstrate that the control time shortens but 
the stability of fidelity weakens with the decrease of the 
antedated time ta- 

Now let us concentrate on the profiles of ilj^ and A 
because they determine the cost of the control and phys- 
ical realizability of the parameters 7 and /3. Figure [3 
illustrates flji due to the condition (v) in Eq. (|18p and 
the conditions in Eq. ([^Uj). while Fig. [51 delineates the 
corresponding detunings. The parameters O/? and A in 
Figs. [3 and [51 arc obviously physical. 




0.0 0.2 0.4 0.6 u.8 1.0 



FIG. 7: (Color Online.) and the energy cost correspond 
to each control passage in Fig. [6l 

To obtain the energy cost of each antedated control 
passage in Fig. [6l we need to integrate each f2/j in Fig. [3 
over time from t = up to the antedated time ta < tf, 
when the population inversion is completed. The legend 
in Fig. [3 presents the result, which shows that the faster 
the antedated control passage the less the energy cost. 
This conclusion may be somewhat rush because we do 
not know yet if the energy cost of each antedated control 
passages under investigation is optimized. It turns out 
to be the case as we will show below. 

The integrals shown in the legend of Fig. [3 suggests 
the existence of a certain lower bound in the energy cost. 
Although it turns out to be true, merely decreasing the 
antedated time ta can hardly lead to the bound, contrary 
to a naive expectation. We point out the reason without 
demonstration. By our study, as soon as ta gets less 
than ~ Itjjl.l, the corresponding 7 will not only cross 
zero at ta but also be below — tt at a later time, which 
brings to fJjj and A extra unphysical behaviors, such as 
discontinuity in $7^ and A and divergence in A, which 
cannot be simultaneously remedied by any form of /3. 
Nevertheless, for any fixed i^, the energy cost admits a 
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Fig. [5] as follows. 



0.2 0.4 0.6 O.S 

t 



FIG. 8: (Color online.) Detuning corresponding to the lEC 
passages in Fig. [S] 



lower bound that is dictated by the value of /3(0), as 
pointed out earlier. 



1.31.7 2.3 



U 3.23 




3.23 



>S'(0) (7i/2tf) 



FIG. 9: (Color Online.) Energy cost as a function of /3(0) (in 
units o{n/2tf). Eadi curve corresponds to the control passage 
in the same color in Fig. [6] with the specified ta- Black dots 
on the curves mark the minima of the curves, whose values 
are recorded in Eq. (pTj). 



To illuminate the relation between 0(0) and the energy 
cost of an antedated control passage, we again consider 
the control passages in Fig_. [6l for each of which we sweep 
0(0) in a wide range, say /3(0) G [0.l7r/2i/, 87r/2</], then 
study the energy cost of the control passage as a function 
of/3(0). The result is plotted in Fig. O Seen in this figure, 
each curve of energy cost has a minimum at certain value 
of /3(0), as marked by the black dots therein. Note that 
${tf) is different from /3(0) by only a minus sign, see 
e.g., Eq. (IT^ (v). We collect the minimal energy costs in 



ta=2tf/7.6 



ta=2tf/7 



ta=tf/3 



Vtndt 



= 3.143, /3(0) = 

tiin 

= 3.144, ^(0) = 

n 

= 3.149, /3(0) = 
= 3.230, /3(0) = 



1.377r 

1.6527r 
2i/ ' 
2.3347r 
2i/ ' 
5.2327r 
2tf 



(21) 



According to Fig. [31 the order of the energy costs in the 
legend of Fig. [7] is no more than a coincidence that arises 
for our unintended choice of /3(0) = 7r/2t/ for those cases. 
Despite this, Eq. (PT|) . which is extracted from Fig. [HI ver- 
ifies that the energy cost of an antedated control passage 
reduces with the antedated time to, as long as the right 
/3(0) is imposed to minimize the energy cost of the con- 
trol passage for a given value of ta- Therefore, antedated 
control passages can be optimized by manipulating both 
ta and /3(0). 

Equation (PT|) also points to the existence of a lower 
bound TT of the minimal energy cost of antedated TEC. In 
other words, an lEC is extremely optimized in terms of 
energy cost when the inversely-engineered Rabi frequency 
rifl, happens to be a 7r-pulse, as /q" f2fldi = tt. Never- 
theless, one must diminish ta further, beyond 2tf/7.Q^ to 
achieve this extremal case. However, this will result in 
unphysical A and before the extreme is hit as ex- 
plained before. 

Now that the energy cost of the antedated control pas- 
sages of lEC can be optimized while the detuning fre- 
quency is physical, stability of the fidelity appears to be 
the only constraint apart from the uncertainty principle 
and the Haniiltonian-switching time. The last condition 
depends on how sharply one can terminate the evolution 
of the system. Given the available experimental power 
nowadays, there would not be any problem in switching 
the Hamiltonian in the scale of nanoseconds in NMR, for 
example. 



IV. CONCLUSIONS AND OUTLOOK 

In conclusion, we have shown by taking the two-level 
system as an example, that lEC indeed offers fast control 
passages — in particular the antedated control passages 
— as shortcuts to adiabatic control of quantum systems. 
To design a LRI for a system, appropriate boundary con- 
ditions on the parameters in the LRI can be imposed to 
meet various goals of the control. In the scheme that is 
based on our polynomial ansatz and the family of bound- 
ary conditions of (3 and 7, an lEC passage is optimal 
in terms of speed and energy cost when the inversely- 
engineered Rabi frequency obtained from the designed 
LRI is approximately a 7r-pulsc. We expect to study 
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the optimization of lEC more systematically by means 
of control theory. All our results apply to pure states, 
which are merely a special case of our general analysis. 
Besides, we note that Carlini et al. [20[ studied time- 
optimal quantum evolution and the optimal Hamiltonian 
for a quantum system via variational principle. It would 
be interesting to see if and how their results are related 
to the optimization of lEC. 

Mixed states are the subject of lEC in this paper, it 
would be intriguing to apply lEC to NMR quantum com- 
putation, in which mixed states are used for computation, 
to further shorten the operation time of the one-qubit 
quantum gates implemented by RF-pulses in NMR. 

It is also natural to extend lEC to multi-level sys- 
tems. To this end, one has to first design reasonable and 
most convenient LRIs for these systems, whose forms are 



rather hard to guess for complicated systems, which is 
under our investigation and will be reported elsewhere. 
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